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Abstract. Hard scattering in p-p collisions, discovered at the CERN ISR in 1972 by the 
. method of leading particles, proved that the partons of Deeply Inelastic Scattering strongly 

■ interacted with each other. Further ISR measurements utilizing inclusive single or pairs of 

\ hadrons established that high pr particles are produced from states with two roughly hack- 

le . to-back jets which are the result of scattering of constituents of the nucleons as described 

' by Quantum Chromodynamics (QCD), which was developed during the course of these 

, measurements. These techniques, which are the only practical method to study hard-scattering 

and jet phenomena in Au+Au central collisions at RHIC energies, are reviewed, as an 
k> ' introduction to present RHIC measurements. 

V 

o 

1. Introduction 

In 1998, at the QCD workshop in Paris, Rolf Baier asked me whether jets could be measured 
in Au+Au collisions because he had a prediction of a QCD medium-effect on color-charged 
partons traversing a hot-dense-medium composed of screened color-charges p. I told him j2] 
that there was a general consensus j3] that for Au+Au central collisions at ^snn = 200 GeV, 
leading particles are the only way to find jets, because in one unit of the nominal jet-finding 
cone, Ar = V / (A?7) 2 + (A0) 2 , there is an estimated tt x ^^fp ~ 375 GeV of energy !(!) The 
good news was that hard-scattering in p-p collisions was originally observed by the method of 
leading particles. 

The other good news was that the PHENIX detector had been designed to make such 
measurements and could identify and separate direct sing le 7 and tt° out to p T > 30 GeV/c 0j. 
It is ironic that the identification of tt° and the separation from direct single 7 out to such a 
large pt in PHENIX was primarily driven by the desire to measure the polarized gluon structure 
function in p-p collisions in the range 0.10 < xt < 0.30 via the longitudinal two-spin asymmetry 
of direct photon production 2 . This illustrates that a good probe of QCD in a fundamental 
system such as p-p collisions, provides a well calibrated probe of QCD in more complicated 
collisions such as Au+Au, should the need arise, which it fortunately did. In the 1980's when 
RHIC was proposed, hard processes were not expected to play a major role. However, in 1998 |2j, 
inspired by Rolf and collaborators, and before them by the work of Gyulassy 5^ and Wang (Hj, 
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I indicated that my best bet on discovering the QGP was to utilize semi-inclusive tt or tt 
production in search for "high pt suppression" . 

It is not generally realized that hard-scattering was discovered — in both Deeply Inelastic 
lepton scattering (DIS) at SLAC J\ and in higher particle production in p-p collisions (EllHltlHj 
at the CERN ISR — before the discovery of QCD. During the period 1972-1982, it was proved by 
single inclusive and two-particle correlation measurements that high pt particles are produced 
in p-p collisions from states with two roughly back-to-back jets which are the result of scattering 
of constituents of the nucleons, at a rate much higher than was predicted from coulomb 
scattering [TT]. This indicated that the partons of DIS interacted strongly (^$> EM) with each 
other. The measurements were consistent in detail with QCD, which was developed during 
this period. In that era, scaling laws relating the transverse momentum (p?>) spectra at different 
values of center-of-mass (cm.) energy, y/s, were the key to understanding the underlying physics. 
Absolute cross sections played a minimal role. Measurements of high pt particle production 
near mid-rapidity, at symmetric p-p or A+A colliders, are ideal for such scaling studies — most 
systematic errors from varying the yfs or the species cancel since the detectors are at rest near 
90° in the cm. system of the reaction and are thus insensitive to longitudinal effects. 

2. Systematics of single particle inclusive production in p-p collisions. 

In p-p collisions, the invariant cross section for non identified charge-averaged hadron production 
at 90° in the cm. system as a function of the transverse momentum pt and cm. energy y/s 
has a characteristic shape (Fig. QJ. There is an exponential (e~ 6pT ) at low pr, which depends 
very little on y/s. This is the soft physics region, where the hadrons are fragments of 'beam 
jets'. At higher pt, there is a power-law tail which depends very strongly on y/s. This is 
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Figure 1. Ed?a/dp s vs. p? at mid-rapidity Figure 2. Data from Fig. Q] plotted as 
as a function of in p— p and p— p collisions. v^(GeV) 6 - 3 x Ed 3 a/dp 3 vs. xt = Zpr/ y/s. 



the hard-scattering region, where the hadrons are fragments of the high pt QCD jets from 
constituent-scattering. 

The hard scattering behavior for the reaction p+p — > C+X is easy to understand from general 
principles proposed by Bjorken and collaborators [121 E] an d subsequent authors [ElEj- Using 
the principle of factorization of the reaction into parton distribution functions for the protons, 
fragmentation functions to particle C for the scattered partons and a short-distance parton- 
parton hard scattering cross section, the invariant cross section for the inclusive reaction, where 
particle C has transverse momentum pt near mid-rapidity, is given by the general 'x^-scaling' 
form [T2| > where xt — 2pT / y/s'- 

& = \n^f) = -Lew. (i) 

dp 6 Pt y/s y/s 

The cross section has 2 factors, a function F (G) which 'scales', i.e. depends only on the ratio of 
momenta; and a dimensioned factor, p^f 1 (y/s~ n ), where n depends on the quantum exchanged 
between constituents. For QED or vector-gluon exchange [IT], n = 4, analogous to the 1/q 4 
form of Rutherford Scattering, while n=8 for the the constituent interchange model (CIM) [T5j . 
quark-meson scattering by the exchange of a quark. When QCD is added to the mix 14 , pure 
n — 4 'vector-gluon' scaling breaks down and n varies according to the xt and yfs regions used 
in the comparison, n — > u(xt, y/s) ^4 — 6. 

It is evident from Fig. [l] that hard-scattering, which is a relatively small component of the 
Pt spectrum at yfs ~ 20 GeV, dominates for pt > 2 GeV/c by nearly 2 orders of magnitude 
at RHIC cm. energies compared to the soft physics e~ 6pT extrapolation. The characteristic yfs 
dependence of the high pt tail is simply explained by the xt scaling of the spectrum. However, it 
is worthwhile to note that it took quite some time for xt scaling with the value of n — 5.1 ± 0.4, 
consistent with QCD, to be observed at the CERN-ISR (Hi El EH- This was due to the so- 
called 'intrinsic' transverse momentum of partons, the "^-effect", which causes a transverse 
momentum imbalance of the outgoing parton-pairs from hard-scattering, making the jets not 
exactly back-to-back in azimuth. The "fcj-effect" acts to broaden the pt spectrum, thus spoiling 
the XT-scaling at values of pt < 7.5 GeV/c, at the ISR, and totally confusing the issue at fixed 
target incident energies of 200-400 GeV (yfs ~ 20 GeV) [E| due to the the relatively steep pt 
spectrum (see Fig. [lj, which results in a relatively strong broadening effect. 

3. Theory and Experiment, c. 1972-82 

The distortion of xt scaling from the QCD prediction due to the fc^-effect is illustrated by 
measurements from the CERN-ISR. The CCOR measurement [HI] (Fig. Et) 3 with a larger 
apparatus and much increased integrated luminosity, extended the tt° measurements of the 
original high pt discovery [HI El QHj to much higher pt- In Fig. [SK, the CCOR 7r° data for 3 
values of are plotted vs xt on a log-log scale. n(xT, yfs) is determined for any 2 values of 
yfs by taking the ratio of invariant cross sections at fixed xt, with results shown in Fig. [Sb: 
n(xT,y/~s) clearly varies with both and xt, it is not a constant. For y/s = 53.1 and 62.4 
GeV, n(xT, y/s) varies from ^ 8 at low xt to ~ 5 at high xt- The fit 16 , for 7.5 < pr < 14.0 
GeV/c, 53.1 < yfi < 62 .4 GeV, is Ed 3 a/dp 3 ~ p^- 1±0A (l - x T ) 12 - 1±0 - 6 (including all systematic 
errors). It is interesting to note that the original ISR high pt discovery in the xt, y/s range 
where n = 8, instead of 4, stimulated the development of the constituent interchange model [T3| . 
which seemed to naturally explain the observed value. The best data at FNAL in 1977 [IZj also 
beautifully showed the CIM scaling with n ~ 8 over the range 0.2 < xt < 0.6, for 200, 300 and 
400 GeV incident energies. However, this effect turned out not to be due to CIM, but to the 
broadening by the /cr-effect. 

An important feature of the scaling analysis (Eq. in determining n(xT, y/s) — is that the 
absolute pt scale uncertainty and many efficiency and acceptance errors cancel! The effect of 
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Figure 3. a) (top-left) Log-log plot of CCOR invariant cross section vs xt — 2pT/yfs\ 
b) (bottom-left) u(xt, y/s) derived from the combinations indicated. There is an additional 
common systematic error of ±0.33 in n. c) (top-right) Invariant cross section for tv° inclusive for 
several ISR experiments, compiled by ABCS Collaboration, d) (bottom-right) ti(xt, y/s) from 
ABCS 52.7, 62.4 data only. There is an additional common systematic error of ±0.7 in n. 



the absoulte scale uncertainty, which is the main systematic error in these experiments, can 
be gauged from Fig. [Sb 18 which shows the tt° cross sections from several experiments. The 
absolute cross sections disagree by factors of ~ 3 for different experiments but the values of 
n(x Tl y/s) for the CCOR [Ej (Fig. Efc>) and ABCS [IHj experiment (Fig. Ell) are in excellent 
agreement due to the cancellation of the systematic errors in each experiment. Thus, while the 
individual ISR experiments each provide a data set with common systematic uncertainties which 
cancel in scaling studies, there is no unique absolute cross section tv° measurement from the ISR 
at 62.4 GeV which can be used as a comparison spectrum for RHIC measurements. 



4. Status of theory and experiment, circa 1982 

Hard-scattering was visible both at ISR and FNAL (Fixed Target) energies via inclusive 
single particle production at large pt > 2-3 GeV/c. Scaling and dimensional arguments 
for plotting data revealed the systematics and underlying physics. The theorists had the 
basic underlying physics correct; but many (inconvenient) details remained to be worked out, 
several by experiment. The transverse momentum imbalance of outgoing parton-pairs, the "fcj-- 
effect", was discovered by experiment QSIEOI? and clarified by Feynman and collaborators pi] . 
The first modern QCD calculation and prediction for high pt single particle inclusive cross 
sections, including non-scaling and initial state radiation was done in 1978, by Jeff Owens and 
collaborators [221 under the assumption that high pt particles are produced from states with 
two roughly back-to-back jets which are the result of scattering of constituents of the nucleons 
(partons). The overall p + p hard-scattering cross section in "leading logarithm" pQCD is the 
sum over parton reactions a + b — > c + d (e.g. g + q — > g + q) at parton-parton center-of-mass 
(cm.) energy = ^Jx\X2S. 

dxidx2dcosv* s ~ 2x\X2 

where f a (xi), fb( x 2), are parton distribution functions, the differential probabilities for partons 
a and b to carry momentum fractions x\ and X2 of their respective protons (e.g. u(x2)), and 
where #* is the scattering angle in the parton-parton cm. system. The characteristic subprocess 
angular distributions, H a6 (cos#*), and the coupling constant, a s (Q 2 ) = ^ ln(Q 2 /A 2 ), are 
fundamental predictions of QCD [231 EH- 

However, jets in 4tt calorimeters at ISR energies or lower are invisible below y/I ~ Et < 25 
GeV [25 . Nevertheless, there were many false claims of jet observation in the period 1977- 
1982 which led to skepticism about jets in hadron collisions, particularly in the USA (SHj- 
A 'phase change' in belief-in-jets was produced by one UA2 event at the 1982 ICHEP in 
Paris 27 J. which, together with the first direct measurement of the QCD constituent-scattering 
angular distribution, E a6 (cos0*) (Eq. EI), using two-particle correlations [2E| 5 presented at the 
same meeting (Fig. 0J), gave universal credibility to the pQCD description of high pt hadron 

physics nana EH- 



5. Almost everything you want to know about jets can be found using 2-particle 
correlations. 

The outgoing jet-pairs of hard-scattering obey the kinematics of elastic-scattering (of partons) 
in a parton-parton cm. frame which is longitudinally moving with rapidity y = l/21n(xi/x2) in 
the p-p cm. frame. Hence, the jet pair formed from the scattered partons should be co-planar 
with the beam axis, with equal and opposite transverse momenta, and thus be back-to-back 
in azimuthal projection. It is not necessary to fully reconstruct the jets in order to measure 
their properties. In many cases two-particle correlations are sufficient to measure the desired 
properties, and in some cases, such as the measurement of the net transverse momentum of a 
jet-pair, may be superior, since the issue of the systematic error caused by missing some of the 
particles in the jet is not-relevant. A helpful property in this regard is the "leading-particle 
effect". Due to the steeply falling power-law transverse momentum spectrum of the scattered 
partons, the inclusive single particle (e.g. tt) spectrum from jet fragmentation is dominated by 
fragments with large z, where z = PTir/PT q is the fragmentation variable. The probability for a 
fragment pion, with momentum fraction z, from a parton with pT q — Prjet is: 
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Figure 4. a) (left 3 panels) CCOR measurement (23 EH| of polar angular distributions of 7r° 
pairs with net pt < 1 GeV/c at mid-rapidity in p-p collisions with y/s = 62.4 GeV for 3 different 
values of tttt invariant mass M n7r . b) (rightmost panel) QCD predictions for H a6 (cos#*) for the 
elastic scattering of gg, qg, qq' , qq, and qq with a 5 ((5 2 ) evolution. 



where D%(z) ~ e~ 6z is the fragmentation function. The change of variables, pT q = Pt^/z, 
dpTq/dpT^lz — 1/^5 then gives the joint probability of a fragment 7r, with transverse momentum 
Pt^ and fragmentation fraction z: 



d 2 a 7T (pT 7r ,z) _ A 



dpr^dz p™~ 



x z m - 2 Dl(z) . (4) 



Thus, the effective fragmentation function, given that a fragment (with p^) is detected, is 
weighted upward in z by a factor z m_2 , where m is the simple power fall-off of the jet invariant 
cross section (i.e. not the u(xt, y/s) of Eq. [l] 32 I). As this property, although general, is 
most useful in studying 'unbiased' away jets using biased trigger jets selected by single particle 
triggers, it was given the unfortunate name 'trigger-bias' jSHj. 

Many ISR experiments provided excellent 2-particle correlation measurements |34j . However, 
the CCOR experiment |35j was the first to provide charged particle measurement with full and 
uniform acceptance over the entire azimuth, with pseudorapidity coverage —0.7 < 77 < 0.7, so 
that the jet structure of high pt scattering could be easily seen and measured. In Fig. Ett,b, the 
azimuthal distributions of associated charged particles relative to a tt° trigger with transverse 
momentum pTt > 7 GeV/c are shown for five intervals of associated particle transverse 
momentum pt- In all cases, strong correlation peaks on flat backgrounds are clearly visible, 
indicating the di-jet structure which is contained in an interval Acj) = ±60° about a direction 
towards and opposite the to trigger for all values of associated pt (> 0.3 GeV/c) shown. The 
width of the peaks about the trigger direction (Fig. EH), or opposite to the trigger (Fig. Et>) 
indicates out-of-plane activity from the individual fragments of jets. The trigger bias was 
directly measured from these data by reconstructing the trigger jet from the associated charged 
particles with pt > 0.3 Gev/c, within Acf) = ±60° from the trigger particle, using the algorithm 
PT)et — PTt + 1-5 ^Zpt cos(A0), where the factor 1.5 corrects the measured charged particles for 
missing neutrals. The measurements of (^trig) — (PTt/PTjet) as a function of pTt for 3 values of 
y/s (Fig. EJ) show the property of xt scaling, which was not expected (SHj- 
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Figure 5. a,b) Azimuthal distributions of charged particles of transverse momentum pt, with 
respect to a trigger tv° with pTt > 7 GeV/c, for 5 intervals of pt- a) (left-most panel) for 
Acj) = ±tt/2 rad about the trigger particle, and b) (middle panel) for Acp = ±tt/2 about tt 
radians (i.e. directly opposite in azimuth) to the trigger. The trigger particle is restricted to 
\rj\ < 0.4, while the associated charged particles are in the range \r]\ < 0.7. c) (right panel) xe 
distributions (see text) corresponding to the data of the center panel. 



Following the analysis of previous CERN-ISR experiments [SIltEj, the away jet azimuthal 
angular distributions of Fig. 03, which should be unbiased, were analyzed in terms of the two 
variables: p G ut = prsin(A0), the out-of-plane transverse momentum of a track; and xe, where: 

PT-PTt _ -prcos(A0) _ z ^ JET T 




XE = ~~\ vT~ = ( 5 ) 

\PTtr PTt ^trig 

z trig ~ PTt/PT]et is the fragmentation variable of the trigger jet, and z is the fragmentation vari- 
able of the away jet. Note that xe would equal the fragmenation fraction z of the away jet, 
for z t rig — ► 1, if the trigger and away jets balanced transverse momentum. The xe distribu- 
tions [SHIISHI for the data of Fig. Efc> are shown in Fig. Et and show the expected fragmentation 
behavior, e~ 6z ~ e ~ 6x E{ztrig) m jf the w idth of the away distributions (Fig. Eb) corresponding to 
the out of plane activity were due entirely to jet fragmentation, then (|sin(A0)|) = (|JtJ/pt) 
would decrease in direct proportion to 1/pt, where is the component of jt in the azimuthal 
plane, since the jet fragmentation transverse momentum, j?>, should be independent of pt- This 
is clearly not the case, as originally shown by the CCHK collaboration [Hih which inspired 
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Figure 6. CCOR j2Hl measurement of (ztrig) as a function of pTt (left) and xxt — ^PTt/Vs 
(right). 



Feynman, Field and Fox (FFF) |21j to introduce, fcj-, the transverse momentum of a parton in a 
nucleon. In this formulation, the net transverse momentum of an outgoing parton pair is y/2kT, 
which is composed of two orthogonal components, v^/ct^ out of the scattering plane, which 
makes the jets acoplanar, i.e. not back-to-back in azimuth, and y/2kT x , along the axis of the 
trigger jet, which makes the jets unequal in energy. Originally, fcy was thought of as having an 
'intrinsic' part from confinement, which would be constant as a function of x and Q 2 , and a part 
from NLO hard-gluon emission, which would vary with x and Q 2 , however now it is explained as 
'resummation' to all orders of QCD FFF [2UE0I gave the approximate formula to derive 
kx from the measurement of p out as a function of xe'- 

(|p out |) 2 = 4[2<KI> 2 + (b'Tj) 2 ] + (\3tJ) 2 ■ (6) 

CCOR 01] used this formula to derive (|&tJ) and (|jrj) as a function of pTt and y/s from the 
data of Fig. (see Fig. Ej). This important result shows that (|jrj) is constant, independent 
of pTt and y/s, as expected for fragmentation, but that varies with both pTt and y/s, 

suggestive of a radiative, rather than an intrinsic origin for fcy. 

6. Final comments 

It should be noted that inclusion of fcy was the key element (221 beyond QCD to explain the 
XT-scaling result of the original CCR 8 high pt discovery and the FNAL (fixed target) 
experiments [17 . More recent FNAL fixed target measurements 02] and many theoretical 
works have used kx as an empirical parameter to improve the comparison of measurement to 
NLO QCD. However, it is important to remember, as illustrated above, that kx is not simply 
a parameter, it can be measured. It is also worthwhile to emphasize that the fc^-effect is 
qualitatively different from NLO QCD. The gaussian nature of kx which is distinctly different 
from the NLO power law tail is crisply illustrated Fig. [HI by measurements of the net-transverse 
momentum distribution of "Drell- Yan" di-muons produced in p-p collisions 43] , a process which 
has zero net pt in LO and diverges in NLO. Another ISR observation j2Hj 5 not much emphasized 
there but relevant at RHIC, is that the measured (^trig) is different for single particle inclusive 
triggers and pair triggers (compare Fig. [Hl-left to Fig. EJ). 
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Figure 7. CCOR 01] measurements of (left) (kr y ) (right) as a function of pTt for 3 

values of y/s. The mean absolute values of the components jr y = jr^ are related to \J (jj) with 
the assumption that jt x and jt v are independent gaussians, with equal r.m.s., which combine 
independently to form j\ = j\ + j\ jl^j . 
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